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ieweeN TRODUCE EON 


Bee LURPOSE AND AYPOTHESIS 


The purpose of this study is to test the applicability 
of the Garwood (1977) one-dimensional bulk mixed layer model 
for hindcasting the warm ocean thermal anomalies found in 
mae Anamoly Dynamics Study (ADS) domain. Descriptions of 
several warm anomalies Suitable for this study are found in 
Elsberry (1984). 

Stringer (1983) conducted a Similar study which involved 
cold anomalies in the ADS domain. He found that the Garwood 
model does a good job of hindcasting cold anomalies in the 
fall-winter period but experiences larger errors in the 
Spring-Summer period. Stringer suggests that uncertainties 
in the heat flux correction fields are fairiy large in the 
Spring and summer and limit the ability of the model to 
hindcast anomalous ocean conditions during those seasons. 
Therefore, a second purpose for this study 1s to further 
examine the heat flux correction fields and to uncover 
factors which may affect the akility of the model to hind- 
cast in tne Spring and summer. 

Three warm anomalies have been chosen for this) study 
that snould satisfactorily test tne Garwood model for the 
spring and summer conditions. The first case involves a 
Warm anomaly starting in January 1976 andending in July 
76. The next two anomaiies extend from June 1977 to 
Mecember 1977 and June _19/¢0 tc Novemper 1978. The time 
periods of these anomalies should adegquatelv test all 
aspects oi the Garwood model for the spring and summer 


hindcasts. 


B. o2UDY DESCREPI AGE 


The North Pacific Ocean Experiment (NORPAX) ADS area 


(Figure 1.1) 1s the region examined in this study. The 
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Figure 1.1 ADS REGION. 


region is bounded by 30N-50N and 130W-160E. Analyses of 
Ships Otero Sport U liek expendable bathythermographs (X Bas) 
provide ocean thermal structure fields on space scales of 
thousands of kilometers and time scales of one month (White 
and Bernstein, 1979). The TRANSPAC analyses provide the 
initial data for the ocean prediction model and also serve 
aS validation fields. 

As briefly described in section A, three warm anomaly 


periods are examined in this study. The first set of warm 


10 


anomalies (WA 1/2) SOceUbowcUrieewanter {Janwary 1976) and 
Seneinues througn July 1976. the second set of warm anoma- 
res, (WA 12/137 liejeeere part ofsa large-scale pattern which 
developed during the summer (June 1977 to December 1977). 
The third set from June 1978 to November 1978 contains an 
extreme, near-surface warm anomaly (WA 23) that forms near 
the western boundary of the domain. All three anomaly cases 
provide excellent examples for tne model hindcast studies 


recommended by Stringer (1983). 
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II. DATA SOURCES AND PREPARATION 


A. DADA SOURCES 


Fleet Numerical Oceanography Center (FNOC) atmospheric 
forcing and TRANSPAC monthly temperature analyses over the 
ADS domain provide the data for this study. The resolution 
is two degrees latitude and five degrees longitude over the 
area bounded by 30N-50N and 130h-160E. 

Atmospheric forcing data were obtained from the arcnived 
FNOC Nortnern Hemisphere analyses and Snort—tegn 
predicticns. The archived data inciude wind speed and direc- 
tion at six-hour intervals and solar flux and total suriaee 
flux at 12-hour intervals. The forcing fields are interpo- 
lated to one-hour intervals for use in the Garwood model. 
This interval is reguired to adequately resoive the oceanic 
response to the diurnal heating cycle. A complete descrip- 
tion of the interpolation procedure can be found in 
Gabilacher (197.9 )i2 

Initialization and verification TRANSPAC data sources 
were described by White and Bernstein (19735 
Ships-of-opportunity XBT observations were made along the 
Shipping tracks between Japan and the west coast or the 
United States. A seasonal and "fair weather" bias, due to 
ships avoiding bad weather, is evident in the nortnern 
portion of the domain, where the data coverage is good in 
the summer and generally poor in tne winter. Elsberry (1984) 
also found tnat the northwest and southwest corners were not 
well sampled and these areas will be excluded in this study. 
These areas will be cross-hatched on appropriate figures. 

Objectively analyzed XBT profiles at depths Ox 
0, 307607 30a oreo and 200m are used jE (O}IE model 


eZ 


initialization. These values are vertically interpolated to 
a 5m Spacing between 0 and 200m to produce tne gridded 
monthly oceanic temperature fields. 

Some observational and physical considerations must be 
kept in mind when studying the subsequent analyses. Elsberry 
(1984) points out that regions north of 45N and south of 35N 
are iess reliabie due to seascnal variations in tne data 
coverage. The baseline study conducted by White and 
Bernstein (1979) indicated that the region west of 1/75W was 
dominated by strong mesoscale (300 km) baroclinic eddies or 
waves, whereas east of 175W the larger scale variability 


dominated. 


Be. PREDICTION MODEL 


The model used for this study is the Garwood (1977) 
second-order closure, bulk model. Advection is not included 
in the model. The Garwood modei is initialized with TRANSPAC 
objectively analyzed temperature fields, which are assumed 
to apply on the 15th of the month. Model output consists of 
vertical temperature profile predictions in 5m increments to 
200m over the entire grid every hour. These hourly temfpera- 
ture profiles are averayed over five days centered oon the 
15th of each month. An additional output is the mixed layer 
temperature at the time of Maximum mixed layer depth for 
each day. Model anomaly Fields are formed by subtracting the 
four-year mean (climatological) fields from tne predicted 
fields at 0, 60 and 120n. 


es 


ITI. WARM ANOMALY DESCRIPTION 


at Ret ee te ee ee eee ee se ee 


A. WARH ANOMALY 1/2 


The warm anomalies (WA) described in this section were 
initially designated and descrited py Eisberry (1984) as WA1 
and WA2. This case persisted for approximately six months, 
aS Shown in Figures 3.1 through 3a 

WA1 1s evident to a depth of 60m during the first month 
of the period (Figure 3.1) in the southeast portion of the 
region. The warm anomaly intensifies to maximum amplitude 
(2.02 C at 34N, 150H) at the surface in February (Figure 
Sepa ss it decreaseS in intensity as the period progresses 
(Figure 3.3 and 3.4). In May, WA1 (Figure 3.5) disappears at 
the surface, but it is still evident at 60m. In June(Figure 
3.6), a warm anomaly ayain appears at the surface in rela- 
tively the same position, but it 1s unclear whether this is 
Stili part of WA1 or is an entirely new feature. 

WAZ 1S barely evident at apforoximately 38N,170E during 
Januapyerigure 3. I) It continues to intensify during the 
period. In March (Figure 3.3), it has become a”significam 
Warm anomaly extending to 120m. In April and May (Figures 
3.4 and 3.5), WA2 appears to tove north and elongate while 
losing some of its structure at depth. Toward the end of the 
period (Figure 3.6), WA2 becomes disorganized at the surface 


but maintains some structure at depths of 60 and 120n. 


Be. WARM ANOMALY 12/13/15 ° 


Warm anomalies 12/13/15 were initially designated and 
described by Elsberry (1984). This case continued for 


approximately six months (Figures 3.7 through 3.12). 
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WA 12 develops in July 1977 in the eastern part 
(oN , 1500) of the region and reaches an anomaly mfagnitude of 
#1.7C (Figure 3.7). Development continues in Auyust (Figure 
3.8) and appears to have reached a depth of approximately 
120m, although the anomaly center is shifted slightly east- 
ward at 60 and 120m. During September and October, WA 12 
weakens and disappears (Figures 3.9 and 3.10). 

WA 13 also deveiops in July 1977 in the region around 
ae, 61 70W (Figure 3.7). In Awvgust 1977, its position and 
movement appear to be influenced by a large cold anomaly at 
approximately 40N, 170W. This ccld anomaly seems to cause WA 
13 to merge with another warm anomaly to the southeast of 
the initial position of WA 13 {Figure 3.8). Evidence of the 
merging in September (Figure 3.%) 1s indicated by the large 
broad warm anomaly in the northeast portion of the region. 
The large anomaly diminishes in size in October and 
November, and it is no longer evident in Decenber (Figures 
Peery and 3.11). 

WA 15 is evident in July 1977 (Figure 3.7) along tne 
extreme western boundary at 4ON, 160E. It continues to 
develop and reaches maximum imtensity of +3.603C in August at 
wa, /OE (Figure 3.8). At this time it also appears to nave 
merged with WA 13, which contrikutes to a large warm anomaly 
area along the northern border. The anomalously high temper- 
atures persist through November and disappear in December 


@ergqures 3.11] and 3.12). 


C. WARM ANOMALY 23 


Warm anomaly 23 was initially designated and descrited 
by Elsberry (1984). This case persisted for approximately 


FOUL montaS aS Seen in Figures 3.13 through 3.16. 


Vs 


During June 19/7@fagqurce 3215 a general warm area at 
the surface has developed on the southwestern boundary of 
the region. By July (Figure 3.14), this area shows a strong 
(approximately +2556) Surface anomaly at the boundary 
region. A slight shift northward is evident in August 
(Figure 3.15), and the temperature anomaly decreases Dy 
approximately 1.5C. By September (Figure 3.16),WA 23 Kas 


disappeared. 


16 


mr a — a a a a ee eS ee 


at eaceseeaO eee eee anaes x 
f 
(A) STO a pects | 
Pome : 
Nesce | . 
A5N Pease 2 } 
RSLS \ A 
Se | 
/ oh 


CA] t> 
on - 
= ae 
os 
24 4 
vest 
*) 
iS \ 
rss 
> 
é aus 
a Os 
7 
- £ 
ae eee GO 
/ 
] 
\ 


Vo eet 
OSS 

CEES 
S59 


ROK 





SON + 


(C) 


MM 
Sane 
SS 


SC 


oe 
& 


de 
Cy) 
ead 
COM OCOLE 
o, 
ee 


40N 


SON 





SON BS 
160E 170E 180 170W 160W {SOW 140W 130N 





Figure 3.1 Temperature (C) anomaly WA 1/2 during 
January 1976 at (A) surface, (B) 60m and (C) 120m. 
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Similar to Figure 3.1 except for February 1976. 
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Similar to Figure 3.1 except for March 1976. 
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Similar to Figure 3.1 except for April 1976. 
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Similar to Figure 3.1 except for May 1976. 
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Figure 3.6 Similar to Figure 3.1 except for June 1976. 
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Similar to Figure 3.7 except for August 1977. 
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Similar to Figure 3. 
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Figure 3.14 Similar to Pigure 3.13 except for July 1978. 
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Begure 3.15 Similar to Figure 3.13 except for August 1978. 
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Figure 3.16 Similar to Figure 3.13 except for September 1978. 
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iV. HODEL HemyDCASTS 
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fee LNITIAL HINDCASTS 


Several model integrations were conducted to test the 
hindcast ability of the Garwood modei for the warm anomaly 
cases described earlier. The correction fieids applied to 
the FNOC forcing fields were described in Elsberry et ai 
932). They correction — Ete) ds, which were derived by 
assuming a locai heat balance during a 36-month period, 
consist of six bi-monthly fields which are to be added to 
the FNOC heat fluxes. The corrections are necessary for the 
heat fluxes to be used for one-dimensional ocean prediction. 
These corrections eliminate a bias of excessive upward 


surface heat flux found in the FNOC fields. 


ieee ode! Handcasts of WA 172 


Generally, hindcasts of WA 1/2 were consistent with 
observations until April. Placement of the two major anoma- 
lies was reasonable and the overall structure for January 
through March was weil simulated. This result is Similar to 
miter OL stringer (1983) for coid anomaly hindcasts prior to 
tne spring transition. In April through Juiy (Figures 4.1 
through 4.4), the model hindcasts tended to be too warm. In 
April, awarm anomaly deveioped at approximately 37N, 1705 
With verified vertical structure to 200 m. The ancmaly 
becomes more intense in May with a temperature of 4.5C at 
the western boundary. Between June and July, temperatures 
continued to increase unrealistically, especially in the 
western portion of the region. The anomalies did penetrate 


to 200 a, aS 1S seen by the high subsurface temperatures. 
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Comparisons with Figures 3.4 through 3.7 Show that theremae 
good agreement in April at a depth of 60n, but the placement 
of the features at the surface and 200 m does not agree with 
the observations. Similar comparisons for May and June show 
even less agreement as the model continues to warm exces- 
Sively. 

An additional model hindcast was conducted using Hay 
as the initial month. Again, the hindcast was too warn, and 
the results very similar to the June initialization. To 
determine the effect of the correction to the forcing field, 
two model integraticns were conducted from January through 
July and May through August without application of the 
correction field. Neither of these model hindcasts agreed 
with the observed anomaly fields in that both predicted a 
warm anomaly in the northwest area anda large area of 
anomalously low temperatures in the entire southern region 
from approximately 40N southward. Further studies with nodi- 
fied correction fields in space and time will be described 


im Section € of this@enapter- 


2- Model Hindcasts of WA 12713715 


Model hindcasts of WA 12/13/15 were generally cooler 
than opservations over much of the area with a pronounced 
Warm region in the southwest pertion of the region. Using 
the June analysis to specify the initial conditions of the 
model (Figure 4.5), it is seen that July (Figure 4.6) is 
already much cooler than climatology in the northeast and 
muco warmer in the southwest. Comparison with observaticns 
for July (Figure 3.7) shows very little correlation between 
the twe flots. The August plot (Figure 4.7) maintains the 
same pattern developed in July and does not bear any resemb- 


lance to the observed July condition (Figure 3.8). 
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The model was then tested with initial conditions 
from March and May to determine if a relationship existed 
between initialization period and accuracy of model hind- 
casts. In both cases, the model predicted development of a 
general warm area in June, although by July there were cold 
anomalies in the northeastern portion of the area. By 
August, the size of the warm area was increaSing (Figures 
mero and &.9). This was similar to the model Output with a 
Mane Initialization. Thus, time of initialization does not 
appear to arfect the output as there are anomalousiy low 
temperatures as early as July. ‘The model was also integrated 
Mernout a Correction to the forcing fieid. The output (not 
Shown) for a model hindcast from June inital conditions 
Showed the entire grid area wasS anomalously cool in July, 
and the temperatures continued to decrease as the period 
progressed. The different variations used to hnindcast WA 
few 7is dia not Signiaficantiy improve any of the hindcasts. 
Thus, the early Summer period appears to pose Serious prob- 


lems for reliable model hindcasts. 


PeeOde! Hind casts of WA 23 


Two model integrations for WA 23 were done with an 
initialization time of January and June. Whee a nee t aia 
zation of June (Figures 4.10 through 4.12), the rodel hind- 
Saots produced a large anomalously cCooi area over most or 
the region with a small warm ancmaly on the southern becrder. 
AS the period progressed throuyh August and September, the 
cold anomalies at the surface intensified and developed an 
Maverified cold pattern. Ccmparison with observations 
(Figures 3.13 through 3.16) indicate a lack of agreement in 
feature placement or feature intensity during all months. 
Temperature structure with depth also did not agree 


mavoraply. 
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General results of each nodel hindcast showed that the 
modei anomalies for WA 1/2 were hindcast well until April. 
However, tne model tended to produce an anomalousiy warm 
condition throughout the remainder of the period. Aodei 
anomalies for WA 12/13/15 were hindcast as too cold as early 
as the second month of the model integration. Model anonma- 
lies for WA 23 were also hindcast as too cold throughout the 
period. 

The three cases have at least one common factor in that 
they span the early Summer nonths of June and July. During 
these months, changes in forcing (heat fluxes and wind 
Speeds) are much more difficult to predict and can not be 
corrected by a bi-monthly average correction technique as 
orevicusly done (Stringer, 1983) for cold anomaly studies. 
The variability in the total heat flux can cause the model 
to over-correct or under-correct in the quantity of) hea 
necessary to simulate accurately the observed conditions for 
the early summer. 

The importance oi correctly determining the heat Stone 
and heat input is discussed by Wyrtki and Ulrich (1382). 
They found that over a period of one month the observed 
change in heat storage can be computed to only an accuracy 
of +/- 100 W/m2?, which 1s as large as the net daily heat 
exchange. This uncertainty makes 2t ditiicmeeece account 
for contributions of advection o© mixing over a Monti 
Intervals. If this magnitude cf uncertainty is encountered 
in actual ocean estimates, then it is very possible that 
errors this large may exist in the FNOC short-tern 
predictions used here for the forcing fields. Uncertainty 
in determining the heat storage and heat input can also 
adversely affect tne depth of the mixed layer as predicted 


by the model. For a Shallowing mixed layer depth with no 
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entrainment, the mixed layer depth ,h, is related to the 
oma necdtmnrl @meneby the Felationship: 
h*Ltu*3 /{gQ,/ (fC,)} where 

L=Obukhov lengtn scale 

u*=Friction velocity 

Cp =Specific heat of water 

QO,=Heat flux 

« =Specific volume 

~ =Density 

g=Gravity 


The corresponding temperature change 1S given by 

dT/dt= (Q #q@ /h. 
The heat flux may be expressed as 

QO,=G. (Correct) +Q.o (error) =Q.+Q,=Q- (1+E), where 
Eis an error bias equal to Q¢g/C.. This now gives 
h=FCpu*3 /fogQe (14E)} = be/ (1+E) 
and with the expression for Q- ylelds 
eae (Oca eee eee iy G) =u 3 | 
Get lt eee t/a 

This expression demonstrates that a net fiux error,d , 
contributes aS a Square because it also contributes to an 
incorrect mixed layer depth. 

The heat £lux corrections by Eisberry et ai.g (1982) were 
computed on a bi-monthly basis and were averayed over tnree 
years. This average correction frield does not appear to 
accurately modify the forcing field in the proper amounts 
for the early summer. To study tnis problem further, the 
Sorrection field for WA 1/2 was modified, and a different 
source of heat fiux values was used to improve nodel 


hindcasts. 
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C., WARIY7Z DETAILED SAUD 


1. One-Dimensionai Heat Budget Correction Field 
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Previous case studies of cold anomalies by Stringer 
(1983) and the initial model hindcasts for warm anomalies in 
this study indicated that the modei hindcasts are degraded 
Ey boundary condition (surface heat flux) errors in the late 
spring and early summer. Since the bi-monthly average 
correction field derived by the method of Elsperry et al. 
(1982) dad not correctly modify tne forcing fields fou yaa 
Spring and early Summer cases, anew correction field was 
derived. This correction field exactly accounts fou eee 
difference between the heat content in the water coiumn toa 
depth of 200 m and the integrated surface heat at a given 
geographical location. The correction is formed from the 


fol owl ag. 


Ey 
T=JSQ,(t) dt 
H= eof ET (2) —T (20g) } az 


I is the integrated surface heat flux, H is the heat 
content, Q;is the surface heat flux, t is time, % is the 15th 
of one month and Tzis the 15th cf the following month. This 


yields the unfiltered exact correction (€); 


€ =(H (+) -H (4) } “17 (b t;) hours 


This field was then filtered using the same procedureas in 
the O§Flginal Correction field 

The one-dimensional heat budget correction field was 
used in wodel hindcasts for May through August 1976. These 
hindcasts were compared with monthly temperature change 
fields. ThiS comparison emphasizes modei temperature 
Changes and is not dependent on climatology as in the 


anomaly comparisons akove. 


Bus: 


For June to May model temperature fields (Figure 
4.13) are close to the observed fields (Figure 4.16). iis 
overall structure, the fieid shows some Similarity to the 
observed pattern, especially in the warming of approximately 
2C in the eastern third of the region. 

The July to June model plot (Figure 4.14) Snows two 
distinct warm features at 42N,175t and 32N,165W which 
fewiatemgreatily from the observed fields (Figure 4.17). 
Laryer temperature difference gradients are more evident in 
the model plot. Locations of maximum wWwarminjy and cooling do 
hot correspond to Similar features in the observed field, 
and it appears that the model is generating temperatures too 
high for the entire region during July. 

The August to July (figure 4.15), Model output 
appears to have even less correlation with the observed 
field (Figure 4.18) than during the previous two montns. The 
two general areaS with iarge deviations noted previously 
again appear. The overall pattern of large gradients in 
both warming and cooling constitutes an unverified hindcast. 

To determine what caused tne model to deviate So 
greatly from observations in July and August, observed and 
model temperature profiles at 42N,175E were analysed 
(Figures 4.19 and 4.20 ). It is seen from the June and July 
observed temperature profiles that a large amount of heat 1s 
added in the water column duriny this two-month period. 
Calculations of the neat storage show that approximately 
3.0x10% cal/fcm? is added, which is abnormaliy large when 
Compared to surrounding values. Analysis of the surface heat 
flux field shows a value of 1.0x10% cal/cm2? which is in yood 
agreement with surrounding points.- Therefore, it appears 
that the large heat content influenced the correction and 
caused the model to over-correct by adding too much heat. 
The reason for such a large neat content could be horizontai 


advection of warmer water, downward verticai velocities, or 
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possible opdservational error. Whatever caused this heat 
excess is not handled well by the model, and must be better 
understood to be corrected. Use oft the ore-dimensional heat 
budget correction field is not surficient for the model to 


produce acceptable hindcasts tor ali points. 


2. Correction to the One-Dimensional Heat Budget 
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Correction Field 


In an attempt to correct the large deviations seen 
in the model using the exact correction field, a correction 
method applied to a specific area was devised. The poi 
(42N,175E) was chosen for this test since there is a large 
temperature difference between July and June and tne nodel 
temperature profile is greatly exaggerated ror July. “aie 
point correction, Ep,was derived by removiny the amount of 
heat inthe mixed layer during the two-month period and 
uSing the residual heat content as an additional correction 
factor to the correction fieid at that point. For (42N 7 17a) 
the following calculations were used: 

Ep =4H2,, ~AHmu. J =29- 13% eal emesis 


where AHzoo =3x10% cal/cn2 and AH, =1.0x103 
cal/fcm2 and t= 720 hours. 
This correction was then applied to the one- 


dimensional heat budget correction field at that point, and 
the model was integrated again. The resuiting nodei tempera- 
ture profiies (Figure 4.21) are much closer to tne observed 
temperature profiles. The July frofile is within 0.5C at the 
surface and is in close agreement down to 80 m. The August 
model prcfile still shows a large amount of heat in the 
column and requires some additional correction to reproduce 
the observed temperature profile. 

Although this procedure did improve the July temper- 
ature profile, it wouid not be useful to apply a similar 


correction at every ;oint without a better understanding of 


4Q 


moe CcduSe Gz the Error {1.e., meglect of advection or okser- 
wa bLOnal GLrOr) . A better nethod for deriving c¢itner the 
correction field or the initial forcing field is needed to 


make better use of tne model aS a prediction tool. 


3. Average Climatological Surface Heat Flux 
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The solar flux values used in the above model hind- 
casts were derived from the FNOC short-term predictions. 
The above attemvts to correct for uncertainties in this 
meela. dalid using Variations of exact corrections did not 
produce overall favorable hinacasts. A climatological solar 
flux estimate will now be used to determine if tsodel 
performance could be improved. 

Climatological heat flux values developed by Wyrtki 
96 7) provide an alternate source for montualy-averaged 
Surface heat flux values in the North Pacific. Values for 
the months orf May through September were extracted from the 
avaliable graphs and then interpolated in time as required 
py the model. 

Monthly climatological heat flux and the average 
heat flux fields from FNOC are compared in Figures 4.22 
merough 4.29. It iS apparent that the FNOC heat fluxes 
contain laryer gradients and more zonal structure. The FNOC 
values tend to be larger when ccmpared to climatology. This 
large deviation from climatology was unexpected and it may 
help explain the difficulty in deriving an objective correc- 
tion to this field, aS waS seen in the earlier model hind- 
Gast S . 

Three model integrations were mace using this clima- 
memwecy. Case i used the totai heat flux field, and it 
neglected tne absorption of shortwave radiation below the 
Surface. That is, all solar radiation was assumed to be 


absorbed at the surface. tigi Case Lisas the solar radiation 
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was added as a separate field that was also derived from the 
Climatological charts of Wyrtki. Thus, below-Suriace 
absorption was considered in Case II. Case III was similar 
to Case II, except it included the diurnal cycle for the 
Gadikation- 

Case I: Plots of monthiy temperature changes for 
June to May (Figure 4.30) were much smoother than the 
previous results uSing the FNOC keat fluxes (Figure 4. 13) 
and temperature magnitudes were much closer to the observed 
values (Figure 4.16). The region has a general warming 
pattern increasing from west to east and north to south with 
a small temperature change maxigum in the southeast corner. 
The July to June plot (Figure 4.31) Shows a larger gradient 
from north to south with a much more intense warming area 
Situated at approximately 32N,170W. The pattern tends to be 
warmer than observations (Figure 4.1/7) with the laryest 
temperature changes shifted slightly south and east of a 
Similar but lower temperature change maximum in the observa- 
tion.The August to July plot (Figure 4.32) also has a 
smoother field and shows a slight decrease in temperature 
Change. The buli'ts-eye pattern apparent in the previous plot 
has been modified and appears as a broad warm area in the 
soutnwestern part of the grid. These temperatures are 
Slightly lower than cbservations (Fijure 4.18). General 
warming throughout the entire region is still evident. 

Although placement of features is not completely 
consistent with observations, the general overall warming 
pattern is verified. Temperature change magnitudes are 
close to observed values. Temperature patterns produced py 
the model tend to be somewhat zcnal due to the zonal pattern 
of the climatological neat flux fields. 

Casewer: Inclusion of a separate solar radiation 
field did modify the temperatures conpared to Case I. The 


June to May plot (Figure4.33) 1s very similar to that in 
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Case I with a reduction of the temperature magnitudes by 
apOEOximately 1.50C. Ther weecomoune | plot (Figure 4.34) 
shows a Similar trend, with temperatures reduced again by 
approximately 1.0C. The August to July plot (Figure 4.35) 
also shows a pattern Similar to that seen in Case [I with 
slightly smaller temperature changes. 

Overall, the inclusion cf a Separate solar radiation 
field produced patterns very Similar to that in Case I. The 
gradients of temperature difference were much smaller, Jonel= 
the observed general warming trend 1S correctly Simulated. 

Caco! (eee UE iNcmevdidatwon to tne Solar radia- 
tion field further modified the temperature difference 
pattern seen in Case Il. Ali three plots (Figures 4.36 
through 4. 38) are very Similar in structure to the Case II 
plots previously discussed. The Juiy to June plot shows the 
greatest modification with an additional decrease of 0.5C in 
the area of the buil's-eye descriped earlier. The general 
warming pattern is stili apparent and, even with the modifi- 
cation of temperatures and weaker gradients, a pattern 
Similar to that of the observations is predicted. 

Climatological heat flux fields do aid the model in 
producing patterns which are clcser to observations than the 
previous attempts using FNOC heat flux fields. HOWeVEL, 
placement of features seen in the obServations is still not 
Hoemdled weli by the use of climatology for the heat flux 


hOorcing. 
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Figure 4.1 
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Model temperature anomaly MA 1/2 during April 1976 


at (A) surface, 


(B) 60m, and (C) 120m 
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Figure 4.2 
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Similar to Figure 4.1 except for May i976. 


45 





SON : : : : 
L60E L7OE 180 {70W 160W 150W 140W 130W 


mm mm i er am a ea mmm i re eee eS 


Figure 4.3 Similar to Figure 4.1 except for June 1976. 
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Figure 4.5 Model temperature anomaly MA 12-15 during 
June 1977 at (A) surface, (B) 60m, and (C) 120m. 
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Figure 4.6 Similar to Figure 4.5 except for July 1977. 
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Figure 4.7 Similar to Figure 4.5 except for August 1977. 
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Figure 4.8 


model initialization in May 1977. 
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Similar to Figure 4.5 except for 





Figure 4.9 


Q 


6D 
SSS 


KxeX™ 


ae 
cote SY 
ee of 
WS 
Yo 


\/ 


Similar to Pigure 


eZ 
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Figure 4.10 Model temperature anomaly MA 23 during July 1978 
at (A) surface, (B) 60m, and (C) 120n. 
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Figure 4.11 Similar to Figure 4.10 except for August 1978. 
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Figure 4.12 Similar to Figure 4.10 except for September 1978 
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Figure 4.13 Model monthly temperature difference MA 1/2 
for June to May with the one-dimensional 
heat budget correction field. 
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Figure 4.14 Similar to Figure 4.13 except for July to June. 
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Figure 4.15 Similar to Figure 4.13 except for August to July. 
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Figure 4.16 Observed monthly temperature difference 


WA 1/2 for June to Hay. 
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Figure 4.17 Similar to Figure 4.16 except for July to June. 
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Figure 4.18 Similar to Figure 4.16 except for August to July. 
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Figure 4.19 Model tal pet tee BOS for May through August 
1976 at oid Deewele 
one-dimensional feat’ kudget correction. 
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Figure 4.20 Similar to Figure 4.19 except for observed values. 
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Figure 4.21 Similar to Figure 4.19 except with added corrections. 
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Figure 4.22 Climatological total heat flux for May. 
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Figure 4.23 FHOC average total heat flux for May. 


65 





x 


a | 






LAC? 


SON oS LLL 


160E 170E 180 170WN 160W 1SQW 140W 130W 


Figure 4.24 Similar to Figure 4.22 except for June. 
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Figure 4.25 Similar to Figure 4.23 except for June. 
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Figure 4.26 Similar to Figure 4.22 except for July. 
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Figure 4.27 Similar to Figure 4.23 except for July. 
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Figure 4.28 Similar to Figure 4.22 except for August. 
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Pigure 4.29 Similar to Figure 4.23 except for August. 
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Figure 4.30 Model monthly temperature difference MA 1/2 
for June to May with climatological heat flux field (Case I). 


69 


ww’ 
< 
OA 
. 
wv 
A 
me 
Y 
«K 
» Y¥ 
s a 
Js ¢ 
= 
« 
rn » 
% 
™ 
¥ 
x. 
x 
4 
¥ 
~~ 
% 
ny; 
x ¢ 
me J 





180 


Similar to Figure 4.30 except for July to June. 


Figure 4.31 
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magure 4.32 Similar to Figure 4.30 except for August to July. 
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Figure 42343 
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climatological radiation field (Case [I). 
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Similar to Figure 4.30 except with a 
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Figure 4.34 Similar to Figure 4.33 except for July to June. 
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Figure 4.35 Similar to Figure 4.33 except for August to July. 
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Figure 4.36 Similar to Figure 4.30 except with a 
diurnal radiation field(Case III). 
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Figure 4.37 Similar to Figure 4.36 except for July to June. 
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Figure 4.38 Similar to Figure 4.36 except for August to July. 
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V. SUMMARY AND CONCLUSIONS 

None of the three warm anomaly (KA) periods chosen for 
this study were hindcast well Dy the sarwood one-dimensional 
model using neat fluxes derived from FNOC. Since all three 
hindcasts seemed to develop difficulties during the early 
Summer, one case (WA 1/2) in 1976 was selected for further 
study. 

A detailed examination of the FNOC surface heat flux 
fields used in previous studies of cold anomalies (Stringer, 
1983) revealed that the early Summer can not be corrected by 
a seasonally varying heat flux correction. Corrections to 
the average heat flux fields did not prove as successful for 
warm ancmalies as for the earlier cold anomaly studies. A 
locai heat budget correction field did not significance 
improve nindcasts. Sohe improvement in specific precbdlen 
areaS waS seen when selective corrections were donc yeas 
this procedure would be of little value in real-time 
predictions. For research purposes, however, selective 
correction is a useful method to diagnose potential model 
difficuities in the hindcast mode of operation. 

Use of climatological heat flux fields in conjunction 
With FNOC winds produced a better result than did the use of 
FNOC's neat fluxes. FNOC heat fluxes and winds and _ the 
various forms used in this study produced general patterns 
which more closely resembled observations. Whenever tempera- 
ture anomalies are primarily attributable to anomalous wind- 
mixing, the use of climatological heat fluxes with actual 
winds (synoptically varying) lay provide reasonable hind- 
casts of thermal structure. However, this climatology, as 
used in this study, still does not account [for the mixed 


layer dynamics attributable to synoptic variabiiity in the 
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heat fluxes. HOrEZCntal abd vertical advection are aiso 
neglected, and may account for some of the differences 
between the observations and the model simulations. Also, 
tne one month resolution afforded py the observations is not 
adequate to verify a model which is most sensitive to 
Synoptic variabiiity in the atmospheric forcing. 

Tiomorldweiaomsioun thatedt iS possible to hindcast 
warm temperature anomalies for the mid-latitudes. The inpor- 
tance of uSing a reasonably correct heat flux field was 
Shown dramatically. This study indicates a need for much 
additional research to improve the surface heat fluxes. 
weeyeeeven wail the model be able to predict (hindcast) 
specific features that are related to the surface thermody- 
hamic boundary conditions. Future three-dimensionai 
modelling will include the advective processes neylected 
here, but the largest difficulty and potential source of 
error still appears to be the specification of the surface 
meat Llux. 
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